Critical exponent for the quantum Hall transition 



Keith Slevin 

Department of Physics, Graduate School of Science, Osaka University, 
1-1 Machikaneyama, Toyonaka, Osaka 560-0043, Japan 



o 
o 

(N 



Tomi Ohtsuki 

Department of Physics, Sophia University, Kioi-cho 7-1, Chiyoda-ku, Tokyo 102-8554, Japan 

(Dated: July 27, 2009) 

We report an estimate v = 2.593 [2.587, 2.598] of the critical exponent of the Chalker-Coddington 
model of the integer quantum Hall effect that is significantly larger than previous numerical estimates 
and in disagreement with experiment. This suggests that models of non-interacting electrons cannot 
explain the critical phenomena of the integer quantum Hall effect. 
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The most important characteristics of a continuous 
phase transition are its critical exponents. The critical 
exponents are a quantitative characteristic of the criti- 
cal fixed point. Agreement of theory with experiment 
implies that the critical fixed point studied theoretically 
is the same fixed point observed in the experiment. For 
the integer quantum Hall transition [l| the value of the 
critical exponent v that describes the divergence of the 
localization length £ of the electrons has now been mea- 
sured to an accuracy of a few percent [3, 0|, so com- 
parison with experiment has become a stringent test of 
theory. One of the most surprising aspects of the integer 
quantum Hall transition is the agreement between exper- 
iment and theory [!, 0, [H, 0] concerning the value of the 
critical exponent. The agreement is surprising because 
the Coulomb interaction between the electrons is not in- 
cluded in the theoretical models [j| 0, [l(J. It suggests 
that the Coulomb interaction is not relevant at the inte- 
ger quantum Hall critical point. Yet, this is immediately 
contradicted when we look at the dynamical critical ex- 
ponent z. For models of non-interacting electrons z is 
exactly 2 [ll|, whereas the measured value is w 1 Q. 

We show in this paper that previous theoretical work 
has significantly underestimated the critical exponent v 
for non-interacting electrons and that there is in fact a 
clear disagreement with experiment. This suggests that 
models of non-interacting electrons cannot explain the 
critical phenomena of the integer quantum Hall effect. 

Non-interacting electron models of the integer quan- 
tum Hall transition are also of interest in their own right. 
In particular, though the critical field theory for such 
models is as yet unknown [13 ]. it has been speculated 
that it should have conformal invariance [l3j|. If so, a 
scaling relation follows [141 ] , which we test by comparing 
recent multi-fractal analyses [HI, Ell with our work. 

The integer quantum Hall effect [l[ occurs in two di- 
mensional electron gases that are subject to a large per- 
pendicular magnetic field. The application of the field 
results in the quantization of the kinetic energy of the 
electrons and the formation of Landau levels. Impurity 
scattering causes the Landau levels to broaden into Lan- 
dau bands. The states at the center of the Landau band 



are critical while other states are localized. A quantum 
Hall transition between quantized values of the Hall resis- 
tance occurs whenever the Fermi energy passes through 
the center of a Landau band [(J [13] ■ The critical expo- 
nent v describes the divergence of the localization length 



near the critical energy. Here, x can be any control pa- 
rameter, for example the Fermi energy, that drives the 
two dimensional electron gas though the transition at 
x c . In common with other continuous phase transitions, 
the value of the critical exponent is expected to exhibit 
a high degree of universality [HI . 

We have performed a finite size scaling (FSS) anal- 
ysis of the quantum Hall transition in the Chalker- 
Coddington network model [HI EH- In this model the 
motion of the electron in a random potential and quan- 
tizing magnetic field is replaced by the transmission of an 
electron through a network of links and nodes. The links 
describe electron motion along lines of constant poten- 
tial and the nodes describe the scattering of electrons at 
saddle points of the potential. The Coulomb interaction 
between the electrons is ignored. 

We use the transfer matrix method to estimate the 
smallest positive Lyapunov exponent of a very long net- 
work consisting of L layers. Each layer consists of two 
sub-layers: one a transverse array of N nodes of type A 
and the other a transverse array of N nodes of type B. A 
transfer matrix XJ of size 2N x 27V relates the N right- 
going and N left-going flux amplitudes at the left of the 
layer to the similar quantities at the right of the layer 
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The explicit form of the transfer matrix is 



Ti = BViAUi. 
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The scattering at nodes of type A is described by the 
matrix 
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Here, O2 is the 2x2 zero matrix, t is the probability 
amplitude with which current from the left is transmitted 
to the right, while r is the amplitude with which current 
from the left is reflected back to the left. Nodes of type 
B are obtained from nodes of type A by a rotation of 90° . 
Imposing periodic boundary condition in the transverse 
direction wc then have 
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The transmission and reflection amplitudes are parame- 
terized by the parameter x according to 



t = (exp (2.x) + 1) 
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The parameter x may be interpreted as the energy of the 
electron measured from the Landau band center scaled 
by the Landau band width. The matrix Ui is a diagonal 
matrix with elements 

(Ul)m,n = <5m,r t exp (vpi, m ) ■ 

and similarly Vi is a diagonal matrix with elements 

( V l)m,n = <W ex P (V*,m) • 

Since the distances between the nodes are random, we 
suppose that the phases are independently and uniformly 
distributed on [0, 2n). (These random numbers were gen- 
erated using the routine ran described in Chapter B7 of 
0.) 

The transfer matrix of the network model is equal to 
the product of the transfer matrices of the layers 



T = ]jT l . 
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From this matrix we can define a hermitian matrix Q by 
n = In (T^T) . 



As a consequence of the conservation of flux the eigen- 
values of this matrix occur in pairs of opposites sign 

{+v\, ■ • • , +v N , -v N , • • • , -vi} , v\ > ■ ■ ■ > v N > . 

We are interested in the smallest positive Lyapunov ex- 
ponents which is obtained in the following limiting pro- 
cedure 



7 = lim 

L — >cx 
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The value obtained is independent of the particular se- 
quence of transfer matrices in the transfer matrix prod- 
uct for almost all sequences. Truncation of the transfer 
matrix multiplication at a large but finite L yields an 
estimate of the Lyapunov exponent with a known preci- 
sion. Repeated QR factorizations are needed to avoid a 
loss of precision due to round off error [20j ]. 

To analyze the simulation data we assume that the 
dimensionless quantity 



r = jN. 



obeys an FSS law 



r = F (N a (x - x c )) . 

Here, a is the reciprocal of the critical exponent v. In 
practice, the data deviate from this law. This is taken 
into account by including corrections to scaling that arise 
from non-linearity of the scaling variables and irrelevant 
scaling variables. In previous work 2l| on the 3D An- 
derson model, we have found that the fitting formula 



r = F (N a u ) + F 1 (N a u Q ) N v ui 



(1) 



works well. Here, uq = uq (x) and u\ = u\ (x) are, re- 
spectively, relevant and irrelevant scaling variables, and 
Fq and Fx are corresponding scaling functions. The expo- 
nent y associated with the irrelevant variable is negative, 
y < 0. The relevant variable is zero at the critical point. 
For the Chalker-Coddington model the critical point is 
known to be exactly x c = and does not need to be 
found by fitting the numerical data, so uq (0) = 0. 

When periodic boundary conditions are imposed, 
swapping t and r, and translating the network so that 
the A and B nodes are interchanged yields an identical 
network. As a result of this symmetry T is an even func- 
tion of x. (Note that this holds only for periodic bound- 
ary conditions.) To reflects this, we make Fq, F\ and u\ 
even, and uq odd, functions. The condition on uq follows 
because the relevant variable is zero at the critical point. 

We are able to fit our numerical data to (pj provided 
that data for smaller system sizes (N < 16) are excluded. 
All the functions appearing in {1} are expanded in Taylor 
series. The results of the fit are summarised in Table 
HI The data and the fit are plotted in Figs. Q] and [5] 



Energy x 



Energy x 



FIG. 1: Simulation data and FSS fit for the Chalker- 
Coddington model. The different lines correspond to different 
numbers of nodes N = 16, 24, 32, 48, 64, 96, 128. 



FIG. 2: The same as Figure 2 but with the focus on data near 
the critical point so as to make the existence of corrections to 
scaling due to an irrelevant scaling variable evident. 



estimate 95% confidence interval 
a 0.3857 [0.3849, 0.3866] 
T c 0.780 [0.767, 0.788] 

y -0.17 [-0.21,-0.14] 

TABLE I: The least squares fit of the FSS model to the sim- 
ulation data. The number of data points is Nr> = 217 and 
the number of parameters is Np = 9. The minimum value 
of x 2 = 199.8 and the goodness of fit p = 0.6. The series 
for Fo, Mo, Fl and m were truncated at orders 6, 3, and 2, 
respectively. 



The confidence intervals have been obtained from Monte 
Carlo simulation [22I ]. 

It is clear from Fig. [5] that corrections to scaling are 
of the order of a few percent for the smallest N. The 
precision of our data is 0.03%, which requires L w 10 s ~ 
10 9 transfer matrix multiplications to achieve. Such high 
precision is required because the dependence on energy 
of r is very weak near the critical point for the available 
N and the critical exponent is, in effect, estimated from 
the variation of the curvature with N. 

We have tested the stability of this fit when the data 
are filtered versus the range of T (Table [TTJ) , the range 
of N (Table IIII[) and also checked for stability against 
increases in the the order of the various Taylor expansions 
(Table lIVp . In all cases there is a large overlap of the 
confidence intervals with the fit in Table U 

Our FSS analysis yields the following best fit value and 
95% confidence interval for the critical exponent 



= 2.593 [2.587,2.598]. 



(2) 



This result is consistent with but also a considerable im- 
provement on the original estimate v = 2.5 ± .5 [101 ] of 
Chalker and Coddington. 



T< 2.0 0.3858 [0.3850,0.3867] 0.779 [0.767,0.788] 
T< 1.5 0.3855 [0.3847,0.3864] 0.781 [0.769,0.790] 
T< 1.0 0.3826 [0.3796,0.3858] 0.789 [0.778,0.797] 

TABLE II: The best fit of a and T c with 95% confidence 
intervals. The data have been filtered using the condition on 
T at the left. 



Our estimate is significantly larger than the oft quoted 
result of v = 2.34 ± .04 of Huckestein and Kramer || 
for a random Landau matrix model. This disagreement 
could be taken as evidence against the universality of the 
critical exponent. However, we feel that a more likely 
explanation is that the precision claimed by Huckestein 
and Kramer is too optimistic. This may be clarified in 
future work. 

Our result is also different from the analytical result 
v = 7/3 of Milnikov and Sokolov However, this value 
is not expected to describe the true critical point but an 
intermediate behavior not too close to the critical point. 

Our result also disagrees with v = 2.37 ± 0.02 of Cain 
et al. [7 1 based on a real space renormalization group ap- 
proach 23j . A more recent calculation [3] on a network 
with a triangular lattice gave v « 2.3 ~ 2.76. Thus, 
the precision claimed by Cain et al. may be question- 
able. Moreover, this approach involves an uncontrolled 
approximation and so the discrepancy may not be signif- 



Ar> 24 0.3854 [0.3845,0.3864] 0.782 [0.763,0.794] 
N > 32 0.3849 [0.3837, 0.3861] 0.787 [0.768, 0.800] 

TABLE III: The best fit of a and L c with 95% confidence 
intervals. The data have been filtered using the condition on 
N at the left. 



4 



F 0.3858 [0.3849,0.3866] 0.779 [0.767,0.788] 

u 0.3853 [0.3843,0.3864] 0.780 [0.767,0.789] 

Fi 0.3828 [0.3754,0.3917] 0.780 [0.767,0.789] 

tti 0.3857 [0.3849,0.3865] 0.780 [0.767,0.788] 

TABLE IV: The best fit of a and T c with 95% confidence 
intervals. The order of the expansion of the indicated function 
has been increased by two compared with TABLE [T] 



quantum Hall effect. Further work is needed to come to 
a definite conclusion concerning the conformal invariance 
of the critical theory of non-interacting electron models 
of the integer quantum Hall transition. 
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icant. 

More important than the discrepancy with previous 
theory is the disagreement of our estimate with recent 
experiments. Li et al 0, Q measured v = 2.38 ± .06 in 
experiments on GaAs-AlGaAs heterostructures. In these 
experiments, an exponent k — 1/vz was measured. Ex- 
traction of the critical exponent v requires an indepen- 
dent measurement of the dynamical exponent z. Li et al 
measured z « 1, but seem to have obtained their result 
for v by supposing z = 1. An unambiguous comparison 
with theory will require a more careful consideration of 
the precision of the estimate of z. 

In our opinion, the coincidence of previous theoretical 
and experimental estimates of v is not significant. Non- 
interacting theory predicts z = 2 in clear disagreement 
with experiment. Calculations [1 11 ] within the Hartree- 
Fock approximation suggest that the observed value of 
z w 1 may be explained by including the Coulomb inter- 
action. We speculate that the discrepancy in the value of 
v that we report here may also be explained in this way. 

Assuming that the critical field theory of the integer 
quantum Hall effect has conformal invariance, it has been 
shown that r c and ao (determined in multi-fractal anal- 
ysis [H, Gil) are related by [3] 



T c = 7r (ao - 2) . 



(3) 



Eq. ^ follows from a conformal mapping between the 
strip used in transfer matrix calculations and the 2D 
plane used in multi-fractal analysis. Inserting our result 
for T c into this formula we obtain 

a = 2.248 [2.244,2.251] . 

This value is not consistent with either 2.2617 ± 0.0006 
reported in [lj| or 2.2596 ± 0.0004 in 0. However, we 
cannot conclude that the critical theory does not have 
conformal invariance because the leading correction to 
scaling decays very slowly with N, as is indicated by the 
small value of the irrelevant exponent y. This does not 
affect the estimation of the critical exponent v but it does 
complicate the estimation of T c and, probably, ao- 

In conclusion, we report an estimate of the critical 
exponent v of the integer quantum Hall effect that is 
significantly larger than both previous theoretical esti- 
mates and, more importantly, experimentally measured 
values. We speculate that models of non-interacting elec- 
trons cannot explain the critical phenomena of the integer 
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